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Abstract. Scalar curvature invariants are studied in type N solutions of vacuum 
Einstein's equations with in general non-vanishing cosmological constant A. Zero- 
order invariants which include only the metric and Weyl (Riemann) tensor either 
vanish, or are constants depending on A. Even all higher-order invariants containing 
covariant derivatives of the Weyl (Riemann) tensor are shown to be trivial if a type 
N spacetime admits a non-expanding and non-twisting null geodesic congruence. 

However, in the case of expanding type N spacetimes we discover a non- vanishing 
scalar invariant which is quartic in the second derivatives of the Riemann tensor. 

We use this invariant to demonstrate that both linearized and the third order 
type N twisting solutions recently discussed in literature contain singularities at large 
distances and thus cannot describe radiation fields outside bounded sources. 



PACS numbers: 0420, 0430 
1. Introduction 

The Petrov algebraic classification of the Weyl tensor and the asymptotic forms 
of radiative fields of spatially bounded sources (peeling theorem) demonstrate that 
solutions of type N play a fundamental role in the theory of gravitational radiation. 
(For a recent review of exact approaches to radiative spacetimes, see e.g. and 
references therein.) 

All solutions of the vacuum Einstein equations with in general non-zero 
cosmological constant A, which are of type N with a non-twisting null congruence, 
are known [ic}| , j|, JlJ]. In the twisting case, only one (Hauser's) solution is available 
||l0f . Although far-zone radiation fields of bounded sources are approximately of type 
N, no known exact type Absolution is asymptotically flat. 

Recently, an interesting discussion appeared in literature jL7), JTs| , || in which 
twisting type N vacuum field equations were solved approximately. Stephani jl8| 
argued, within a linear theory, that no solutions regular outside a bounded source, 
exist. Finley, Plebahski and Przanowski H iterated the solution up to the third order 
and concluded that their iterative procedure leads to the regular solutions, so that " it 
seems that the twisting, type N fields can describe a radiation field outside a bounded 
sources." 

In neither of these works, however, the arguments are really compelling since 
singularities are not analyzed properly. (For general treatments of singularities, see 
e.g. Q.) The authors study only the behaviour of invariants with respect to gauge 
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transformations which leave invariant coordinate system and field equations for type 
N metrics because for type N vacuum spacetimes with A all scalar invariants of the 
Weyl or Riemann tensor are trivial - they either vanish or are constants depending on 
A. 

It is known that for a vacuum pp-wave, which is a very special case of a type A?" 
metric, in fact all curvature invariants of any order, i.e., invariants depending also 
on covariant derivatives of the Weyl (or Riemann) tensor, vanish [Q. However, 
nothing appears to be known about higher-order invariants in more general type N 
spacetimes. The problem of such higher-order invariants is addressed in the present 
paperf . 

We define an invariant of the metric gapix 1 ) of the order k as a non-constant 
scalar function I(g a /3, g a /3,-yn ■ ■ ■ , 9a/3,n...~f k+2 ) which satisfies 

I(9aP,9al3,7i> ■ ■ ■ > g a f3,- n ...-y k + 2 ) = ^(<? Q| 3, 3 a/3, 71 1 • • ■ > 9 a/3 ,71. ..7^+2) 

under a spacetime diffeomorphism x a — > x' a — x' a {x^). It can be proved (see e.g. 
H ) that any invariant of the order k depends on the metric and the Riemann tensor 
Raf3-/S and its covariant derivatives of the order < k: 

I = I(9a/3> RaP-ySi ■ ■ ■ j Raf3-y8-ei...£ k ) ■ (1-2) 

It has been known from 1902 what is the maximal number, I(k,n), of functionally 
independent invariants of the order k in a Riemannian space of dimension n. Denoting 



D(k;n)=I(k;n)-I(k-l;n), (1.3) 
Haskins Q found that 

mn) = (3)^ forfc = 0, (1.4) 

Hence, in a 4-dimensional spacetime there exist in general 14 functionally independent 
invariants of the order zero which depend only on g a p and Rap^s, 60 invariants 
depending on g a p, R a p-yS, Ra(3~/S;p, 126 invariants containing also R a /3-yS;pa-, etc. 
Clearly, when the derivatives of the Riemann tensor are included, the number of 
independent invariants grows rapidly. All 14 invariants of the zero order were explicitly 
given in , their spinor equivalents can be found in |2C| ] . 

In Section 2 we first briefly review basic definitions and relations of the two- 
component spinor formalism and Newman-Penrose formalism flljj , which will be 
needed later. Then we demonstrate the generally known fact that in vacuum type 
N spacetimes with A = all zero-order invariants vanish. If A 7^ 0, some invariants 
are non- vanishing but they are just constants depending on A. Using spinor formalism 
we prove a helpful Lemma 1 on the properties of the invariants constructed from the 
derivatives of the Weyl tensor. 

In Section 3 we specialize to the type A vacuum solutions with A, admitting a non- 
vanishing and non-twisting null geodesic congruence (called Kundt's class in @). We 
prove that in these spacetimes all invariants constructed from the Weyl tensor and its 

f It is partially based on the thesis |15| 
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covariant derivatives of arbitrary order vanish. The invariants constructed from just 
the Riemann tensor are all constants depending on A. 

Expanding and non-twisting type N solutions are discussed in the first part of 
Section 4. Using again extensively the spinor and Newman-Penrose formalism we find 
that all the first-order invariants vanish. However, the formalism indicates that a non- 
vanishing second-order invariant may exist. A number of attempts have eventually 
led to the non-vanishing invariant 

I = R al3 ^ 5 ^ 'R a ^^- £l j > R XfJ ' pl/ '' J7 R\p p S;oT- (1.6) 

The zero and first-order invariants vanish also in the twisting case. The invariant ( |l.(i| ) 
remains non- vanishing. 

This invariant is then, in Section 5, used to analyze the nature of approximate 
solutions |Q , H we mentioned earlier. We find that Stephani's conclusion, based 
on the linearized theory, remains true for the third order solution obtained by Finlcy 
et al : both solutions contain singularities at large r. This raises more doubts about 
the physical meaning of type N twisting solutions as describing radiation fields outside 
bounded sources; nevertheless a definitive statement can be made only when a general, 
exact solution is found. 

The technique by which we arrive at the non-vanishing invariant of the second 
order could probably be employed also in more complicated cases of algebraically 
special spacetimes. Scalar invariants obtained might play a role not only in classical 
relativity but also in a quantum context. 



2. Higher-order curvature invariants: general properties 

In this chapter we prove a lemma about the properties of the curvature invariants 
of higher order in Petrov type N spacetimes. The lemma will be used in the next 
section in which curvature invariants in the specific classes of type N spacetimes will 
be analyzed. First, we have to summarize some basic notations and relations which 
will be needed later. 



2.1. Spinors and the Newman-Penrose formalism 

Spinors and the Newman-Penrose formalism have been reviewed by many authors (see 
e.g. fl(i[| , [pj| , fig}]). Here we give only a very brief summary. 

Consider a null congruence of geodesies with the tangent null vector l a , 

i a r = o, (2.i) 

which is affinely parametrized, 

rf. a = 0. (2.2) 



(Hereafter, we assume Eq. ( j2.2[) to be satisfied.) The congruence is characterized by 
expansion 9, twist ui and shear \o~\ given by 

=\l°\ a , ( 2 -3) 



-IfrflW ~ 6 2 • (2-5) 
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In algebraically special spacetimes the shear vanishes. A 2-component spinor field 
o A € W (A = 1, 2, W is a 2-dimensional complex vector space) can be associated with 
a null vector field l a in a standard way, 

l a < — ► o A d A , (2.6) 

where o A € W', and another spinor field, i A , exists such that together with o A it forms 
a spinor basis satisfying 

o A u A = e AB o a l b = 1, o A o A = l a l a = 0, (2.7) 

where the Levi-Civita alternating symbol e AB plays the role of the metric in spinor 
calculus. 

The complex null tetrad, {l a , n a , m a , fh a }, in which l a is introduced by fl2.6|) and 



n" < — > t"L", m" < — ► o"t", m a < — > l a o a , (2-8) 
satisfies the usual relations 

l a n a = -fh a m a = 1. (2.9) 
The (covariant) derivative operator V Q can be expressed in the form 

V Q = n a D + l a A — rh a S — m Q 5, (2-10) 

where 

D = l a V a , A = n a V a , S = m a V a , S = m a V a . (2.11) 

In terms of the covariant derivative with spinor indices, 

= u% V Q (2.12) 

(a a nd a % are proportional to the unit and Pauli matrices), we have equivalently to 
Eqs. (EOT)) 

D = o A d A V AA , A = tVv^, 
S = o A l A V M , 6 = t VV^. 



The Newman-Penrose twelve complex scalar spin coefficients k, e, 7r . . . are defined as 
frame components of the covariant derivatives of the null-tetrad vectors. Our notation 
follows that customarily used in literature jL2], [^9|. Some details are spelled out in 
Appendix A. 

The Weyl tensor Ca^s in the spinor form is given by 

C a (3jS i > ^jbcd £ab£cd ~\~ ^ Abcd £ A b£co 5 (2-14) 

where ^ AB cr> — ^(abcd) ■ In the Newman-Penrose formalism the Weyl tensor 
is described by five complex scalar quantities ^o, ^l, ^4 given by the 

projections of ^abcd onto basis spinors o a ,l a . The Riemann tensor R a p-ys in the 
spinor form is given by 

Raf3yS * y X ABCD E Ab Ecd + X AB( jo £ ab £cd ,^ -^gs 

+ 3 > ABCD &AB&CD ~\~ ^ABCD &AB&CD1 
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where 



R 
12 



-AD C BC ) ) 



R is the scalar curvature, and the spinor & AB cn — &(ab)(c£>) — ®abcd corresponds 
to the traceless Ricci tensor S a p = Ra/3 — jRga/3'- 



'S a b 



(2.16) 



In vacuum spacetimes with generally non-vanishing cosmological constant A, 

R = 4A, (2.17) 

and 



0. 



(2.18) 



Bianchi identities connect the derivatives of \l/'s with the \& 's themselves and with 
the Newman-Penrose spin coefficients. All Bianchi identities, the relations giving the 
Riemann (Weyl) tensor in te rms o f the spin coefficients and the commutation relations 
for the derivative operators ( 2.11 ) are explicitly written down in |l(J, for example. 

Here we are interested in the vacuum spacetimes of Petrov type N with in ge nera l 
A ^ 0. Let o A be the 4-fold principal null spinor of ^abcd, and l a satisfies (2.7). 
Then we have 



*o = *i = *2 = *3 = 0, 
^ABCD = ^ao a o b o c o d , 



(2.19) 
(2.20) 



and R and $ AB cb are given by Eqs. ( 2.17 ) and ( 2.18 ). The Bianchi identities reduce 
to 



D* 4 = (p - 4e)* 4 , 
= (r - 4/3)^4, 
a = k = 0. 



(2.21) 
(2.22) 
(2.23) 



Although under the conditions (2.17)^(2.20) and (2.23) the remaining Newman- 
Penrose equations take much simpler form than in general case, they still represent 
the set of 21 equations. These are given in the Appendix. 

As is well-known, transformations preserving the direction of l a can be divided 
into two subclasses: 



(i) Null rotations around vector l a 



r = i c 

I A , 



cl a 



n' a =n a + cm a + cm a + ccl c 



o 



i' A = L A +CO A . 

(ii) Boosts in the (l a ,n a ) - plane and spacelike rotations in the (m Q ,m Q ) - plane 



(2.24) 



l' a = a 2 l a 



n = a n 



(2.25) 



where z 



ae 



() 



In the following we shall in particular need the behaviour of various quantities under 



the pure constant boosts in the (l a , n a ) - plane (a =const, 9 = in (2.25)). 
If a quantity fl transforms under the boosts as 

ft' = a q n, (2.26) 

the number q is called the boost-weight of fi. We write b(fl) to denote the boost- 



weight. First notice that b{o A ) — 1, b(t 



-1, and recall that operators D,A,S and 



S are also weighted. The boost-weights of the Newman-Penrose coefficients (NP) and 
operators (OP) are summarized in the following Table 1: 



x e NP 


xeOP 


b(x) = b(x) 


K 




4 


spa 


D 


2 


a (3 7T T 


5 





7 A H 


A 


-2 


V 




-4 



We shall also need the behaviour of ^4: 
Higher-order curvature invariants 



(2.27) 



As noted in the Introduction, in a general spacetime there exist 14 independent 
invariants of the zero order, i.e., invariants depending only on the metric and the 
Riemann tensor. In terms of the spinors determining the Riemann tensor according 
to relations ( [2.14 ), these invariants can be constructed as products of the form 

ifi if^ABMN \r -yABMN tf^GDAB if, \rRSKL if, /p, M N 

^ABMN * , A-ABMN A ) * 1?RSAB A ^KLMN ^ GD , 



etc. Since, however, in the vacuum Petrov type TV spacetimes with A equations (2.17)- 
(2.20) are satisfied, it is easy to see that among the zero-order invariants 9 vanish and 
5 are dependent just on the value of A (as, e.g., ( [2.17 )). Hence, we have to turn to 
the invariants of higher order. 

First, let us decompose the spinor derivative V "*" . . .\7 ClXl (^>40 A o B o c o D ) into 
the spinor basis of the appropriate spinor space. Let us write 
yy[p,fc] 



W x W x 
p 



x W x W x W x 
k 



x W. 



Thus, VVtP' k l is 2 p+k - dimensional complex vector space whose basis spinors 



B 



[P,k] j^[pM 



B^p+l can be constructed from the tensorial products of o A and l a 
of the form f^ 1 £ A 2 . . . £p p Xf 1 Aj 2 ■ ■ • A^ fc , where £ A< 's are o Ai or i Ai and \* 3 's are 
o x j or I x j . Then the decomposition of the n-th spinor derivative of ^> i o A o B o c o D reads 



yc n x n ..V c A(* 4 o A B C O i3 ) = ^ Ci-B. 



[n+4,n] 



(2.28) 



We shall need to know some restrictions on the coefficients Cj rather than their 
specific values. Let us first note that the coefficients Cj must all be the sums of the 
terms of the form 



A1A2 . . . X n ^4 



(2.29) 
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where Xi e NP or S OP. This can easily be seen by regarding the well-known 
decompositions of V** in terms of the NP operators (2.11) and of V^o 3 , V 4 *^ 8 in 
terms of the basis spinors o, t and the NP coefficients. 

It is now useful to introduce a simple notation. If a product of basis spinors has 
the form Y = p Al . .. o A ™\o^ . o x ^yi , . ^ , then P (Y) = mi +m 2 

mi ma Tii ri2 

will denote the number of o A 's and o A 's , and Pi(Y) — n\ + n 2 of i A, s and l A, s which 
are contained in Y. 

Now it is easy to see (essentially as a consequence of (2.7)), that if a spinor 
S A 1 ... Am x 1 ...x k e >v[™,fc] has the form 



,k] 



(2.30) 



then all invari ants formed from the products of S " must vanish provided that the 
coefficients in ( 2.30 ) are such that Sj = for all % for which 



Po(B. 



m,k] \ 



< P L (B\ m ' k] ). 



(2.31) 



This observation enables us to prove the following 
Lemma 1: 



Let an invariant constructed from the products of the spinors 

yc„x„ _ _ _ yCiXj (\j/ 40 A s c D ) j for fixed n, be non-vanishing. Then there exists a non- 
vanishing quantity X\X 2 ■ ■ ■ X n 1 i l 4, Xi e NP U OP, such that 

n n 

b(X x X 2 . . . X n * 4 ) - Y, b ( X ^ + 6 (*4) > 0, i.e., J2 ^ 4 - 



Proof: The n-th spinor derivative i s of the form ( 2.28 ), where the coefficients Cj 
are sums of the terms of the form ( 2.29 ). According to the observation above, an 
invariant formed out of these derivatives will be non-vanishing only if there exists 
Q 7^ such that P (B\ n+4 ' n ^) < P t (_B,[ n+4,n '). Since the basis spinors have the boost- 
weight d = b(B\ n+4M ) = P (B\ n+4 ' n] )-P L (B\ n+4 ' n] ), and since from dJD and Q 
it follows that 6(cj) = —d, we see that Cj must satisfy the condition 6(cj) > 0. Hence, 
there must exist a non-vanishing quantity X\X 2 ■ ■ ■ X„^ 4) Xi £ NP U OP, such that 
b{X x X 2 . . . X n f 4 ) > 0. 



3. Non-expanding and non-twisting solutions 

Type N vacuum solutions with A admitting a non-expanding and non-twisting null 
geodesic congruence belong to Kundt's class (see |lo| , ch. 27, for details). Since 
9 = lu = (cf. Eqs. (2.3), (2.4)), the NP coefficient 

p = 9 + iuj = 0. (3.1) 

In this chapter we shall prove that in this class all the curvature invariants of any 
order vanish (or are constants determined by A), generalizing so results of || and Q 
where only the plane- wave metrics are considered. 

Choose the null tetrad parallelly propagated along the null congruence determined 
by the multiple principal null direction of the Weyl tensor and parametrized by an 
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affine parameter. Thus, only *f> 4 ^ (cf. equation (2.19)) and NP coefficients 
(7 = k = 7t = £ = 0. The NP equations, given in Appendix for a general type N 
vacuum spacetime with A, then simplify considerably. From all the NP equations we 
shall only need those containing the operator D: 



Dt 


= 0, 


Da 


= 0, 


D(3 


-o, 


L>7 


= ra + f/3 


DX 


= 0, 




R 


Dfi 


~ 12' 


Dv 


= fjJL-\- t\, 



El 

24' 



(3.2) 



and the commutators 

(AD - DA) = (7 + 7)£> - t6 - f 6, (3.3) 
(SD-DS) = (a + /?)£>. (3.4) 

We shall also need the Bianchi identity (2.21) which now simply reduces to 

£>*4 = 0. (3.5) 

The simple form of the above equations suggests the following notation: let T k be 
the set of functions / such that / £ T k -^==> D k f = 0. From the NP equations (3.2) 
we easily find that 

a, (3, t, X G J 7 !, 
7, fj, e T 2 , v e T%. 

Using this, and employing the equations (3.3) and (3.4) for the commutators, we can 
prove the following 

Lemma 2: 

Let / e T k . Then (i) Sf £ T k , Sf £ T k , 
(ii) Af£T k+1 . 

Proof: (i) can easily be proven by induction. Applying (3.4) on /i, we immediately 
get D5fi = => Sfi £ T\. Assuming then Sf k £ T k and applying (3.4) on f k +i, we 
find D k+1 8f k+1 = => 8f k+1 £ T k+1 . 

In order to prove (ii), we first show, by using Leibniz's formula, that D k+1 (f2f k ) = 
for all k > 1. Then (ii) can be proven by induction similarly as in (i) (the commutator 
(3.3) now being used instead of (3.4)). 

It will now be useful to associate the number p with any NP coefficient X which 
will indicate the behaviour of X under the action of the operator D. Let 

p{X) = k-l, if X £ T k but X <£F k -i. (3.7) 

If X £ OP, i.e., X is one of the NP operators, we define (being motivated by Lemma 2) 
p(A) = 1, p(d) = 0, p(D) = -1. 

The values of p for all relevant quantities are summarized in the following Table 2: 
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X e NP 


X eOP 


P(X) 


V 




2 


1 V 


A 


1 


a (3 t X 


5 







D 


-1 



The indicators p enable us to formulate easily Lemma 3: 

Consider (as in Lemma 1) a quantity X1X2 ■ ■ ■ X n where X{ € NP U OP. 

If E7=iP( X i) < 0, then X X X 2 . . . X n ^ 4 = 0. 

Proof: From the Bianchi identity ( |3.5[ ) we have — 0. Regarding Table 2 we 

observe that the condition X)"=i P(-^i) < requires that with any Xi having a positive 
p{X{), the operator D must appear at least p(A^)-times among X\ . . . X n since only 
p(D) is negative (for example, with any v, D must appear at least two-times ); an 
additional D has then to enter Xi . . . X n in order that 53i=iP(^j) < 0- This results 
in Xi . . . X n $ 4 = 0. 

Combining Lemma 1 and Lemma 3 we can now prove the following 

Proposition 1: 

In type N vacuum spacetimes with A admitting a non-expanding and non-twisting 
null geodesic congruence all n-th order invariants formed from the products of spinors 
V 2 ™*" . . . V ClXl {^iO A o B o c o D ), with n arbitrary but fixed, vanish. 
Proof: According to Lemmas 1 and 3, a non- vanishing invariant will exist only if there 
are X; e NP U OP, i = l...n, such that 



Y.KXi) >4 and ^ppQ)>0. 



(3.8) 



The values of these sums depend on how many times the specific NP coefficient or 
NP operator enter Xi . . . X n . Let mi denote the number of the coefficients v and 
v, rri2 — of 7,7, jjL, and /2, 7713 — of A and A, 7714 — of a, a, (3,(3, r and f, k\ — the 
number of operators D, k 2 — of S and 5, and £3 — of A, which enter X\ . . . X n . Then, 
regarding Tables 1 and 2, we find out that the inequalities (^^) read 

— 4toi — 2m2 — 2to3 — 2&3 + 2ki > 4, and 2mi + 7712 + k$ — k\ > 0. 

Combining the last two inequalities, we obtain 7773 < —2. This is impossible since all 
777,'s and /c's must be non-negative. 

From Lemmas 1 and 3, and from the last proof it follows that all the coefficients 
c-i in the decomposition (2.28) which multiply basis spinors satisfying the condition 
(2.31), must necessarily vanish. Hence, all non-vanishing terms in the decomposition 
(2.28), with 77 arbitrary, contain a larger number of spinors o A than of l a . Therefore, 



regarding (2.7) it is evident that even all invariants formed from the products of the 
derivatives (|2.28) with different 77's must necessarily vanish. Recalling the relations 



(2.14) and (2.20) between the Weyl tensor C a /3 7 5, the Weyl spinor ^abcd and the 
NP scalar ^4, we can now formulate our two basic propositions. 

Proposition 2: 

In type N vacuum spacetimes with A, admitting a non-expanding and non-twisting 
null geodesic congruence, all invariants constructed from the Weyl tensor and its co- 
variant derivatives of arbitrary order vanish. 
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Next, recall that in vacuum spacetimes with A Einstein's equations imply R a p = 
A<7 a/ g , so that 

C a 0-/S;s = Ra/3fS;E- (3-9) 

Then, regarding the spinor form ( [2.15 ) of the Riemann tensor we see that all invariants 



constructed from the Riemann tensor and its covariant derivatives of arbitrary order 
can be expressed in terms of the spinors 

^ 4 o A o B o c o D + ^(e AC e BD + e AD e BC ) 7 (3.10) 
V 3 ^(9 4 o A o B o a o D ), (3.11) 



V°"*" . . . V lXl {^ i0 A o B o c o D ). (3.12) 

Applying the considerations above and recalling that e AB = o a l b — o b l a , we easily 
make sure that all terms in the invariants containing ^! i o A o B o c o D or their derivatives 
of arbitrary order, vanish. The only non-vanishing quantities can be formed from the 
constant term j(e AC e BD + e AD e BC ) and are dependent on A only. We thus finally 
arrive at 
Proposition 3: 

In type N vacuum spacetimes with A, admitting a non-expanding and non-twisting 
null geodesic congruence, all invariants constructed from the products of the Riemann 
tensor and its covariant derivatives of arbitrary order vanish provided they contain 
a derivative of the Riemann tensor. The invariants constructed from the Riemann 
tensor itself are all constants depending on A. 

If A = 0, the Riemann tensor is equal to the Weyl tensor and all the invariants vanish 
by Proposition 2. 



4. Expanding and twisting solutions 

This section is divided into two parts. In the first, we analyze type A spacetimes with 
A with expanding but non-twisting null congruences. We shall show that invariants 
of the zero and first order vanish (or are constants determined by A) as in the non- 
expanding and non-twisting case. However, we will succeed in finding a non- vanishing 
invariant of the second order depending on the expansion and on ^4. In the second 
part, we shall show that the invariants of the zero and first order again vanish, and 
we shall demonstrate how the non-vanishing invariant is modified when there is a 
non-zero twist. 

Before we turn to the details we wish to point out the main reason why, with 
the expansion present, a non- vanishing invariant may exist: it is due to the Newman- 
Penrose equation for the expansion, Dp = p 2 . Now D n p ^ for any n, and one even 
cannot formulate Lemma 3, for example. 



4--1. Non-twisting case 

All metrics of vacuum type N spacetimes with A were given by Garci'a-Diaz and 
Plebahski ||. In their coordinates, suitable for our purposes, the metric reads 

ds 2 = 2v 2 d£,di + 2vAd£du + 2vAd£du + 2ipdudv + 2{AA + TpB)du 2 , (4.1) 
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where 



A = e£-vf, 



• 4 = —dv — Bdu. 



B= -e + |(/, t +/, f ) + ^ 2 V, 

^ = 1 + e££ , (4.2) 
e = - 1, 0, +1. 
It is useful to choose the null tetrad corresponding to the forms 
uj 1 = vdt; + Adu, 
uJ 2 = vd£ + Adu, 
uj 3, = 4>du, 

UJ 

Using this tetrad one obtains the NP coefficients as follows: 

a = K = £ = ir = \ = a = f3 = 0, 

f>t l AV ^ 1 (A A\ 

7 = ^ + T' T = -~v> p = v> (4 - 4) 

2i/> 6 ' 2ipv 

In the following calculations the specific forms of most of the NP coefficients are 
not important - we need to know only the value of p = 1/v, and the fact which of the 
coefficients vanishes. Then we can write down all the NP equations by specializing 
the NP equations in the Appendix to the present case. The Bianchi identities (2.21), 
(2.22) simplify now to the form 

£1*4 = ,0*4, <5*4 = r* 4 . (4.5) 

Using these and the NP equations we can write the first spinor derivative of the Weyl 
spinor as 

V EF {^ 4 o A o B o c o D ) = - {A + A- f )^ i S ABCDE o F 

(4.6) 

+ (54-4 S ABCDE l F + $ 4 (™ F - P l F )S ABCDE , 

[0] [1] 

where we have introduced quantities S ABCDE which denote the symmetrized products 

of the basis spinors o A and i A , and the subscript " gives the number of i A 's, entering 

S ". (For example, S AB — o a l b +o b l a .) As usually, S " denotes the complex conjugate 

[i] 



to S'", with [j] being the number of l A, s. From equation (4.6) it is then clear that all 
the first-order invariants of the Weyl tensor vanish since there are not enough t's to 
be combined with o's, as was discussed in detail in the previous section. 

Now consider the second derivatives. Since calculations become rather lengthy we 
turn to the computer algebra package Maple V . Again using the NP equations and 
Bianchi identities, we arrive at the following results: 

V H V EF (^iO A o B o c o D ) = AS ABCDEG S FH + BS ABCDEG S FH 

[0] [0] [1] [0] 

+ C s ABCDEG S FH + V s ABCDEG S FH + 8 S ABCDEG S FH (a 7) 

[0] [1] [2] [0] [1] [1] V • ) 

_|_ QABCDEG CjFH _|_ Q gABCOEG _|_ QjABC JDEG 

[2] [1] [1] [2] [2] [2] 
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where functions A, B,C ... read 

A = |(A + 97 + 7)A + 4(57 + 7 + A)7 - 5ti/- 1/6^4 
<5{(A-47)* 4 }, 
-j{(A + 4 7 )* 4 }, 



B 
C 



V = 2* 4 t 2 , (4.8) 
£ = (DA + 4p 7 )*4, 
T = -2* 4 pr, 

g = -2^*4, 

H = 2* 4 p 2 . 



The above expression (4.7) indicates that a non-vanishing invariant may exist. 
Looking at the last term in equation (4.7), which is proportional to the function TL 1 
we see that it contains, in contrast to all other terms, the same number of o's and 
t's. Consequently, a combination of such terms may give a non-vanishing result. One 
can make sure, however, that simple squares or even cubes of such terms do not work. 
After a number of unsuccessful attempts to construct a non-vanishing expression we 
arrived at the invariant 



or, regarding (3.9), equivalently 

j = C^^C^.^C^P^Cxpp^r. (4.10) 
In terms of spinors employed to arrive at the invariant, we obtain 



/ = W 2 H 2 S BGDH S BGDH S B OOH S BGDH , (4.11) 
[2] [2] [2] [2] 

where we have used the relation S ABCDEF i E i F — S ABCD . Since a straightforward 

[2] [2] 

calculation gives S BGDH S BGDH = 6, we find 

[2] [2] 

/ = IUH 2 H 2 , (4.12) 
or, regarding (4.8) and (4.4), we finally obtain 

/ = 9(2p) 8 * 4 2 *4, (4.13) 

and 

^=144^&. (4.14) 
4-2. Twisting case 

Twisting type N vacuum spacetimes with A are not known, except for the Hauser 
solution. With non- vanishing twist, p becomes complex, p = 8 + icu, and in contrast to 
non-twisting spacetimes, the NP coefficients a, (3, A in general do not vanish. Hence, 
the NP equations become much more complicated. However, inspecting the relation 
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generalizing ( |4.6| ), we again find out that the zero and first order invariants vanish. 
The calculations of the second derivatives (4.7) are only feasible using a computer 
algebra package (here again Maple V ). Functions A, B,C . . . in (4.8) become much 
more lengthy, however, the most relevant function, Ti, remains the same. Therefore, 
using (4.12) we find that the invariant (4.9), resp. ( 4.10| ), becomes 



9(2p) 4 (2p) 4 * 4 2 * 



(4.15) 



Before moving on to possible applications of the invariants, let us note that 
their form in the NP formalism (not in terms of the Weyl or Riemann tensor) 
could have been anticipated by considering the behaviour of the NP quantities under 
the transformations (2.24) and ( 2.25| ). Taking into account only undifferentiated 
NP quantities and assuming a type N spacetime with p / 0, one finds that the 
quantities which are invariant under ( [2.24 ) and (2.25) must be functionally dependent 
on p 2 p 2 ^^l^. All such invariants can thus be written in terms of the invariant 



(4.9). Of course, expression (4.9) is invariant under any tetrad and coordinate 
transformation. 



5. Applications 



The invariant (|4.9| ) can be used to study the occurrence of singularities in various 
expanding type N spacetimes. 



5.1. Expanding, non-twisting spacetimes 



For example, in the Garcia Diaz - Plebahski spacetimes with the metric (4.1), we 



find that the invariant (4.9), which now becomes equal to the expression (4.14), is 
diverging if for f^fggg ^ there is (i) v — 0, (ii) ip = 1 + e££ = 0, i.e., if e = — 1 
and = 1. It also diverges when v and ijj are finite but /|^^/|^ diverges. If / is 
a polynomial quadratic in £, the spacetime can be shown to be flat or de Sitter (if 
A ^ 0). If, however, 

/(£, u) = co (it) + ci(u)£ + . . . + c n {u)C, 
where n > 6, then / is singular at £ = co, £ = oo for each e = —1, 0, 1. 



5.2. Type N twisting spacetimes 

Recently two papers appeared discussing the physical meaning of the vacuum Einstein 
equations of Petrov type N with expanding and twisting null congruence. Stephani 
|l8| has argued, within the framework of the linearized theory, that these solutions 
contain singular lines ('pipes') in space which at any time extend arbitrarily far away 
from a possible insular source. A few months ago, Finley, Plebahski and Przanowski 
H , using an iterative approach, tried to resolve " Stephani's paradox" by constructing 
the solution up to the third order in their approximation scheme. They conclude that 
"up to the third order, there do exist acceptable, regular solutions". What does our 
invariant ( 4.15j ) say about the nature of these approximate solutions? 

The standard form of the type N twisting vacuum solutions (see e.g. fyl) reads 
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(using here the usual signature +2) 
ds 2 = 2uj 1 uj 2 - 2w 3 w 4 , 

U) = = LU , 

P'P (5.1) 

lu 3 = du + Ld( + Ld(, 

lu 4 = dr + Wd( + WdC 

where the real function P(C, C u ) an d complex function L(£,£,u) appearing in the 
metric satisfy the relations 

p=-^-, 2 i £ = P 2 («9L-dL), 

t (5-2) 

W = — + idZ , d = d c - Ld u . 
P 

The function p is indeed the NP coefficient used before. Introducing a real function 
V by putting P = V u , the field equations can be written in the form 

(ddddV) tU = p- 1 (ddV), u (BBv), u , 

lm{ddddV) = , d^- 1 {BBV) <u ] = 0. 
The coordinate system and field equations are invariant under gauge transformations 
@, @ 

C' = /(C), u' = F((X,u), r' = rF-\ (5.4) 
The NP component of the Weyl tensor is given by 

* 4 = P 2 P d u [p- x {BBv), u ]. (5.5) 

Under the gauge transformation (^J), the quantities appearing above transform as 
follows: 

p' = F, u p, E = 

d' = r x d, p'=F; u x \r\p, 

i (5-6) 
L = f (Lf, u — -F,c)) 

Using the expression ( p 7 ^ ) and ( |5.5| ) for p and we can easily calculate the 
invariant ( |f .15| ) to obtain 

I = 9(2p) 4 (2p)% 4 2 § 2 = 2304 (d u {p-\ddV), u ]y \du\p- 1 {ddV) , u ]f '.(5.7) 

It can easily be checked, by using relations ( |5.6| ) , that ([T?]) does not change under the 
gauge transformation (5.4). Of course, in its original forms ([1.9[), (4.10), / is invariant 



under any coordinate transformation. 

Now Stephani jl8| found the general solution of the linearized field equations in 
the form (using Stephani's notation) 

P = l + f, 

T = B(C n i c(u ' c " } i C " 2J( "' C) fog i p (5 ' 8) 

{ ' (l + CC/2) 2 2(l + CC/2) 2 (l + CC/2) 
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where functions C(u, £) and D(u, £) are arbitrary, -B(£, £) has to satisfy Im(J5 ^) = 0. 
In order that spacetime is not flat, i.e. \&4 7^ 0, the condition 

D, uuCCC + (5.9) 

has to be satisfied. 



Using Stephani's solution (^8) in (5/7), we find the invariant to be 



1 = ™^ D X<CuuD^ uu . (5.10) 



Defining C(u, () = C(u, £)£ 2 and D(u, £) = D(u,()C 2 , Stephani calculates the 
expression 



P 2 L uc - = 2CC ^; C) + ?C U , - 2 <°^ - < 2 D uC , (5.11) 

• uC 1 + CC/2 ' c 1 + CC/2 



which is invariant under the gauge transformation ( p.4| ). Since functions C and D are 
analytic in £ and £, respectively, they will have singularities in the plane (C>C)j i- c - 
on the sphere (recall the standard convention £ = \/2tan |e* v ). The expression (5.11) 
will be regular only if C jU ~ £ _1 and D u ~ £ _1 . However, in this case the spacetime 
is flat. 

Clearly, function D(u,£) has singularities in the plane (£, C) and the invariant 
( 5.10 ) will diverge unless D tUU ~ a£ 2 + 6£ + c. (Notice that if D ~ £ 3 than the 
invariant diverges due to P 8 .) However, this implies flat spacetime again. 

As mentioned above, Finley, Plebahski and Przanowski [gj constructed the twisting 
and diverging type N solution up to the third order of an iteration procedure. Their 
final expressions for functions entering our invariant (5.7) read 



w /i + iccYd 3 / 



(5.12) 



r V 



C / du 3 ' 

where aW is a complex constant, = f^ 2 \u) + /( 3 )(u), f^(u), f^(u) being 

(21 

arbitrary functions of it. The curvature is non-vanishing so far as /.„„„ 7^ 0. Both 
quantities, E/r and if/r 2 (for if see equa tion (5.9) in Q), considered in 0, which 
are in vari ant under gauge transformations (|5.4[), a re indeed regular; the NP coefficient 
p (cf. (t).2j)) entering our curvature invariant ( 1.15| ) is also regular. However, (ty^^) 2 
is not regular on the (C,C) ~~ sphere unless a^fl ««« = which corresponds to flat 
spacetime. Indeed, (^^i) 2 ~ (1 + ^CC) 8 / (CC) 4 diverges at CC ~~ * 00 (which in 
standard convention, with £ = \/2e lip tan |, corresponds to = 7r). 

Therefore, we find that Stephani's conclusion based on the linearized theory 
remains true for the third order solution analyzed by Finley, Plebahski and 
Przanowski. Although this raises more doubts about the interpretation of type N 
twisting solutions as representing radiation fields outside bounded sources, solutions 
of full Einstein's equations may perhaps bring us surprises. 
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Appendix A. The Newman Penrose equations in type N spacetimes 

The NP coefficients are defined by the following table (see e.g. [^9|), in which "V" 
denotes respectively D, A, 6, S; the first line gives the definition in terms of the basis 
spinors, the second in terms of the null tetrad. 



V 


o A Vo A 


o A Vi A = l a Vo a 


L A Vi A 




m a Vl a 


\{n a Vl a - m a Vm a ) 


-m°Vn a 


D 


K 


e 


7T 


A 


T 


7 


V 


S 


a 





/" 


5 


P 


a 


A 



In the vacuum type N spacetimes the null tetrad can be chosen so that 

*o = *i = *2 = *3 = 0, 
®aAbb =0, a = k = 0, 

the NP equations are 

p 2 + (e + e)p, 
(r + n)p + (e - s)t, 
[p + e- 2e)a - f3e + (e + p)tt, 
(p-e)P- (a-7f)e, 



Dp 
Dt 

Da — Se 
D/3-Se 

Dj - Ae 

DX-Sn 

Dp — Sir 

Dv - Att 

AX-lv 

S P 

Sa -8/3 

SX — Sp 
Sv — Ap 
5-/-A/3 
St 

Ap - St 
Aa — ^7 



R 

(t + 7f)a + (f + 7r)/9 — (e + e)7 - (7 + 7)e + T7r - — , 
pX + ir 2 + (a- (3)ir - (3e - e)A, 



pp 



(e + e)p - 7r(a — /3) + 



R 
12' 



(tt + f )p + (7f + t)A + (7 - 7)7r - (3e + e)i/, 
- (p + p)A - (37 - 7)A + (3a + (3 + tt - f)i> - * 4 , 
p(a + /3) + (p - p)r, 

fip + aa + f3/3 — 2a[3 + 7(p — p) + — A) + 



(A.l) 



24' 



(p - p> + (p - p)tt + p(a + (3) + X(a - 3/3), 
(p 2 + AA) + (7 + 7)p - Ptt + (r - 3/3 - a>, 
(r — a — /3)7 + pr — — /3(j — 7 — p) + aA, 
Ap + (t + /3 - a)r, 

R 

— pp+ {(3 — a — t)t + (7 + 7)p - — , 

(p + e)v - (r + /3)A + (7 - p)a + (0 - f) 7 , 
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and commutators 



(AD - DA) 
(5D - D5) 
(5 A - AS) 

(55 - 55) 



(7 + 7)£> + (e + e)A - (r + n)5 
(a + /3-n)D-(p + e-e)5, 



— vD + (t — a — (3) A + X5+ (p, 
(p-p)D + (p-p)A-(a-(3)5 



(f + n)5, 



7 + 7)5, 
{fi - a)5. 



(A.2) 



References 

[1] Bicak J 1997 Radiative spacetimes: exact approaches, in General Relativity and Gravitational 
Radiation, Proc. of the Les Houches School 1995, eds. Marck J A and Lasota J P, Cambridge 
University Press. 

[2] Clarke C J S 1993 The Analysis of Space-Time Singularities, Cambridge University Press, 
Cambridge. 

[3] Garcia— Diaz A and Plcbanski J F 1981 All nontwisting N's with cosmological constant, J. Math. 
Phys. 22, 2655. 

[4] Ellis G F R and Schmidt B G 1977 Singular Space-Times, Gen. Rel. and Grav. 8, 915. 

[5] Finlcy J D, Plebanski J F and Przanowski M 1997 An iterative approach to twisting and 

diverging, type-N, vacuum Einstein equations: a (third-order) resolution of Stcphani's 

"paradox", Class. Quantum Grav. 14, 487. 
[6] Gchcniau J and Debever R 1956 Les invariants dc courbure de l'espace de Riemann a quatrc 

dimensions, Acad. Roy. Belg. Bull. CI. Sci. 42, 114. 
[7] Haskins C N 1902 On the Invariants of Quadratic Differential Forms, Trans. Amer. Math. Soc. 

3, 71. 

[8] Jordan P, Ehlcrs J and Kundt W 1960, Abh. Akad. Wiss. Mainz, Math. Nat, No 2, p. 77 
[9] Klingcr R 1989 Differentialinvarianten und Symmetrieklassen von Tensoren, Technischcn 
Hochschule Zrich, Zrich. 

[10] Kramer D, Stephani H, MacCallum M and Herlt E 1980 Exact solutions of Einstein's field 

equations, Cambridge University Press, Cambridge. 
[11] Newman E T and Penrose R 1962 An approach to gravitational radiation by a method of spin 

coefficients, J. Math. Phys. 3, 566. 
[12] Penrose R and Rindler W 1984, 1986 Spinors and space-time 1: Two spinor calculus and 

relativistic fields, 2: Spinor & Twistor methods in space-time geometry, Cambridge University 

Press, Cambridge. 

[13] Pirani F A E 1965 Introduction to gravitational radiation theory, in Brandeis (1964) Lectures 

on General Relativity, vol. 1, Prentice-Hall, Englewood Cliffs., New Jersey. 
[14] Podolsky J 1993 On exact radiative space-times with cosmological constant, PhD dissertation, 

Department of Theoretical Physics, Charles University, Prague. 
[15] Pravda V 1998 PhD dissertation, Department of Theoretical Physics, Charles University, Prague. 
[16] Schmidt H J 1994 Why do all the curvature invariants of a gravitational wave vanish ?, in New 

Frontiers in Gravitation, ed. G. A. Sardanashvili, Hadronic Press, Palm Harbor, pp. 337-344. 
[17] Stephani H 1980 Some new algebraically special radiation field solution connected with Hauser's 

type N vacuum solution, Abstracts of GR9, Jena, F. Schiller University, Jena. 
[18] Stephani H 1993 A note on the solution of the diverging, twisting type N vacuum field equations, 

Class. Quantum Grav. 10, 2187. 
[19] Stewart J 1990 Advanced general relativity, Cambridge University Press, Cambridge. 
[20] Witten L 1959 Invariants of General Relativity and the classification of spaces, Phys. Rev. 113, 



357. 



